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Abstract 

When the singlet-doublet fermion dark matter model is extended with additional Z 2 -odd 
real singlet scalars, neutrino masses and mixings can be generated at one-loop level. In this 
work, we discuss the salient features arising from the combination of the two resulting simplified 
dark matter models. When the Z 2 -lightest odd particle is a scalar singlet, Br(^ —)• ej) could be 
measurable provided that the singlet-doublet fermion mixing is small enough. In this scenario, 
also the new decay channels of vector-like fermions into scalars can generate interesting leptonic 
plus missing transverse energy signals at the LHC. On the other hand, in the case of doublet-like 
fermion dark matter, scalar coannihilations lead to an increase in the relic density which allow 
to lower the bound of doublet-like fermion dark matter. 


1 Introduction 

In view of the lack of signals of new physics in strong production at the LHC, there is a growing 
interest in simplified models where the production of new particles is only through electroweak 
processes, with lesser constraints from LHC limits. In particular, there are simple standard model 
(SM) extensions with dark matter (DM) candidates, such as the singlet scalar dark matter (SSDM) 
model [1, 2, 3], or the singlet-doublet fermion dark matter (SDFDM) model [4, 5, 6, 7, 8, 9]. In this 
kind of models, the prospects for signals at LHC are in general limited because of the softness of 
final SM particles coming from the small charged to neutral mass gaps of the new particles, which 
is usually required to obtain the proper relic density. In this sense, the addition of new particles, 
motivated for example by neutrino physics, could open new detection possibilities, either trough new 
decay channels or additional mixings which increase the mass gaps. 
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On those lines, scotogenic models [10], featnring nentrino masses snppressed by the same mecha¬ 
nism that stabilizes dark matter, are being thoroughly studied with specihc predictions in almost all 
the current terrestrial and satellite detector experiments (For a review see for example [11]). The sim¬ 
plest models correspond to extensions of the inert doublet model [12, 13] with extra singlet or triplet 
fermions. Recently, the full list of 35 scotogenic models with neutrino masses at one-loop [14, 15]^, 
and at most triplet representations of SU{2)l, was presented in [17] (and partially in [18]). The next 
to simplest scotogenic model is possibly the one where the role of the singlet fermions is played by 
singlet scalars, and the role of the scalar inert doublet is played by a vector-like doublet fermion. 
One additional singlet fermion is required to generate neutrino masses at one-loop level. This kind 
of extension of the singlet dark matter model is labeled as the model T13A with a = 0 in [17]. The 
extra fermion, required in order to have radiative neutrino masses, can be the singlet in the SDFDM 
model. 

In the simplest scotogenic model [10], singlet fermion dark matter is possible but quite restricted 
by lepton flavor violation (LFV) [19, 20]. In contrast, we will show that in the present model the 
region of the parameter space, corresponding to fermion dark matter, is well below the present and 
near future constraints on Br(/i —>■ ey). 

On the other hand, when the lightest Z 2 -odd particle (LOP) is one of the scalar singlets, in 
the regions of the parameter space compatible with constraints from LFV, we could have promising 
signals at colliders, thanks to the electroweak production of fermion doublets and possible large 
branchings into charged leptons. 

The dark matter phenomenology of both the SSDM and SDFDM models has been extensively 
studied in the literature and recently revisited in [21]. Here we consider the possible effect of coan¬ 
nihilations with the scalar singlets for fermion dark matter. We will see that these coannihilations 
tend to increase the relic density of dark matter and may modify the viable parameter space of the 
model. Specihcally, they allow to reduce the lower bound on the mass of the doublet-like dark matter 
particle from around 1100 GeV down to about 900 GeV. 

The rest of the paper is organized as follows. In the next section, we present the model. Our main 
results are presented in Sections 3 to 6 where we describe the correlation between the generation of 
neutrino masses and lepton flavor violation, new signals at colliders in the case of scalar dark matter, 
and new coannihilation possibilities in the case of singlet-doublet fermion dark matter. Finally, 
in section 7 we present our conclusions. In the Appendix we present the analytic diagonalization 
formulae for the mass matrix of neutral fermions. 


2 The model 

The particle content of the model consists of two S'f/(2)L-doublets of Weyl fermions Ru, Rd with 
opposite hypercharges; one singlet Weyl fermion N of zero hypercharge, and a set of real scalar 
singlets Sa also of zero hypercharge. All of them are odd under one imposed Z 2 symmetry, under 
which the SM particles are even. The new particle content is summarized in Table 1. The most 

^The general realization of the Weinberg operator at two-loops have been undertaken in [16] 
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Symbol 

{SU{2)L,U{l)y) 

^2 

Spin 


(1.0) 

— 

0 

N 

(1.0) 

— 

1/2 

Rui 

(2.+1/2) 

— 

1/2 

Rd 

(2. -1/2) 

— 

1/2 


Table 1: a-set of scalars and Weyl fermions of the model. 


general Z 2 -invariant Lagrangian is given by 

L =£sm + MoeatRA - “ ^deabH^R^d^ - XueabH^R^N + h.c 

- S^Sp + X^JeabH'^H’^S^Sp + X^^sSaS^S^Ss] , 


( 1 ) 


where Li are the lepton donblets, we have defined the new S't/(2)2,-donblets in terms of left-handed 
Weyl fermions as 


Rd — 




and W = (O {h + v)/V2)'^ as the SM Higgs donblet with H = ia 2 H* and v = 246 GeV. In the scalar 
potential, we assnme that the M| matrix has only positive entries and (Tf|)Q ,_3 + Af^n^ = Ofor a ^ (3, 
which means Sa are mass eigenstates with masses and ms^ < the 


other hand, the Z 2 -odd fermion spectrnm is composed by a charged Dirac fermion = (-0^, 
with a tree level mass m^± = Mo, and three Majorana fermions arisen from the mixtnre between the 
nentral parts of the SU{2)l donblets and the singlet fermion. By defining the fermion basis throngh 

T 


the vector H = , the nentral fermion mass matrix reads 


where 




/ 

j —mx cos /3 
y mx sin /3 


—mx cos 13 
0 

-Mo 


mx sin (3\ 
-Mo , 
0 / 


mx = 


Xv 

71’ 



tan (3 



( 2 ) 

( 3 ) 


The specihc signs on the right hand side of eq. (1) were chosen so that the terms in the mass matrix 
follow the same conventions as in the nentralino mass matrix [22]. From this follows that the 
supersymmetric case corresponds to the pure bino-higgsino limit with mx = mzsmOw leading to 
A = g'/\/2. The Majorana fermion mass eigenstates X = (xi,X 2 ,X 3 )^ cire obtained through the 
rotation matrix N as H = NX such that 


NTJVIxN = (4) 

with = Diag(m]^,m^,mg) and m)j being the corresponding masses (no mass ordering is im¬ 

plied). In which follows we assume CP invariance and therefore N can be chosen real. The analytical 
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Figure 1: One-loop Weyl-spinor Feynman rules [25] for the contributions to a neutrino mass, with 
three Majorana fermions n = 1, 2, 3, and a singlet scalar S. 


diagonalization of the neutral fermion mass matrix is carried out in Appendix A. For the subsequent 
analysis, it will be convenient to have some approximate expressions in the limit of small doublet- 
fermion mixing (m^ Mq, Mjsi). Expanding the analytical expressions for the eigensystem of eq. (4) 
given in Appendix A, up to order the fermion masses are 


rrq =Mm + 


Me sin (2/3) -|- Mat 

M% - Ml 


7ni + O {ml) 


V ,, sin(2/3) -f 1 , / A\ 

mf =Mo + 77TTV-^ ("^a) 


rug = — Md + 


2 {Md — Mjsf) 

sin(2/3) - 1 


ml + O {ml) . 


2 {Mjy -|- Mpf) 

Approximate expressions for the mixing matrix are also given in that Appendix. 


(5) 


3 One-loop neutrino masses 

By assigning a null lepton number to the new helds in the modeF, the only lepton-number violating 
term in the Lagrangian eq. (1) is the one with coupling hia- Hence, the introduction of real singlet 
scalars allows to generate non-zero neutrino masses at one-loop level through the diagram shown in 
hgure 1. The resulting one-loop neutrino mass matrix was presented in the interaction basis in [15] 
and [23], and more recently in the limit = 0 and Mjv —)■ 0 in [24]. Instead, we work out the 
calculations in the more convenient mass-eigenstate basis, in which the neutrino mass matrix takes 
the form 


h. 




Ibvr^ 


m^ 


Bo (O; 


m 


Xn 


, rric 




n=l 


( 6 ) 


where Bq is the Bq Passarino-Veltman function [26] and (Nmn) are matrix elements of 

the rotation matrix N. By using the identity 

3 

E(^3„)"< = (M*)33 = 0, (7) 

n=l 

^If complex singlets instead real singlets are considered, an accidentally conserved lepton number would have been 
obtained in the Lagrangian, and such a case vanishing neutrino masses are expected. 
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we obtain the expected cancellation of divergent terms coming from the mass independent term in 
Bq, leading to the hnite nentrino mass matrix 


f , ( 8 ) 

a n=l 

^ ^ ^ia-^a^ja (9) 

a 

= (hAhT),, , (10) 

with / (mi, m 2 ) = (m^ Inm^ — m^ lnm 2 )/(m^ — m^), A = Diag (Ai, A 2 , A 3 ) and 

1 ^ 

Aa =Y^5^(^3n)^rnx„/(m5.,"ixu)- (H) 

n=l 


The flavor structure of the neutrino mass matrix Mt)-, given by eq. (10), allows us to express the 
Yukawa couplings in terms of the neutrino oscillation observables (ensuring the proper compatibility 
with them) through the Casas-Ibarra parametrization introduced in [27, 28]. Thus, by using an 
arbitrary complex orthogonal rotation matrix TZ, the Yukawa couplings hi^ are given by 

h^ = D^7?,D^17t^ 

where = Diag \/m ^), = Diag ^v^A 7 \ y/A^\ j and U is the PMNS [29] 

neutrino mixing matrix. Henceforth we will consider the case of three scalar singlets, a = 1,2,3, 
where the Yukawa couplings take the form 

^ + ^/m:^na2U*2 + 

^ 

In the above equation, the 3x3 matrix TZ can be casted in terms of three rotation angles 6 ^ 23 , ^ 13 , ^ 12 , 
which are assumed to be real. It is worth mentioning that for the case two scalar singlets a = 1, 2 
a viable scenario is also possible with the remarks that one massless neutrino is obtained. To fully 
exploit the generality of hi^ couplings obtained from ( 12 ), we stick to the case with three scalar 
singlets. 

In summary, the set of input parameters of the model are the scalar masses m^^, Mjv, M^, A, 
tan/3, the lightest neutrino mass mi,i, the three rotation angles present in IZ and Af^^. With no 
lose of generality we assume for the latter to be small Af^ < 0 . 01 , except for the case of scalar dark 
matter where is set to give the proper relic density. 

In order to have an approximate expression for in terms of this set of input parameters, we 
can use the identity (7) to obtain 

Aa [f{ms^,m\) - /(ms,,m^)] + Nl^ml [f{ms^,ml) - /(ms,,m^)]} . 

^The couplings irrelevant for phenomenological purposes. 
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(b) 


SIM, A = 5 X 10-3 



SIM, A = 5 X 10-3 


O a 
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tan j5 


Figure 2: tan/3 dependence of (a) and (b) Br(/i —)■ 07 ), for the set of input masses in eq. (14) 
with A = 5 X 10-3. 


The expression for the matrix elements at O {m\) are given in the Appendix A. Since and 
~ cire already we can use the leading order values for the other 

masses and mixings parameters to obtain 


A ^ 

'■a ~ 


IOtt^ 


[/(ms„, Mtv) - fims^^Mo)] + -Md [f{ms^,m^) - m^)] \ + O [m\) . 


With the last two approximate formulas for masses in (5), and the Al|^ mixing in (26), we have 


1 r- 2 
IGtt — 5 - 

mi 


Md cos f3 + Mat sin f3 


Ml-Ml 


+ 


Mn [/(m5„, Mtv) - /(m^,, Md)] 
Ml [Md sin (2/3) + M^] 


{Ml-Ml) [Ml-miy 


Ml-ml^ log 


Ml 

mi 


+ 1 


+ O (mi) . (13) 


To illustrate the dependence in tan/3 of A^, we consider the following set of input masses (SIM) 
compatible with singlet scalar dark matter: 


ms^ = 60 GeV ms 2 =800 GeV =1 500 GeV 

mAr = 100 GeV m£)=550GeV. (14) 

The results for A = 5 x 10-3 are shown in hgure 2(a). For large values of tan/3, the A^ are positive. 
However, there are specihc values of tan /3 for which each A^ goes to zero and turn to negative values 
as illustrated by the red lines in the plot. The specihc point with /3 = vr/O is illustrated by the yellow 
stars in the hgure. 


4 Lepton flavor violation 

The size of the lepton havor violation (LFV) is controlled by the lepton number violating couplings 
hia- From the approximate expression for Aq, in (13) and the analysis of the previous section, we 
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will show that these couplings are inversely related to the Yukawa coupling strength A. Since in 
SDFDM the observed dark matter abundance is typically obtained for A > 0.1 [9], the lepton flavor 
observables are not expected to give better constraints than the obtained from direct detection 
experiments. Therefore we will focus our discussion of LFV in regions of the parameter space where 
Si is the dark matter candidate. 

It is well known LFV processes put severe constraints on the LFV couplings and in general on 
the model’s parameter space. One of the most restrictive LFV processes is the radiative muon decay 
/i —>■ ey, which in the present model is mediated by same particles present in the internal lines of the 
one-loop neutrino mass diagram. The corresponding expression for the branching ratio reads 


Br(p —>■ ey) 


3 CTem 

4 IQirG'jp 


Y^hi^- 


{Ml 


ML) 


h. 


m 


2a 




(15) 


where 


F{x) 


— 6x‘^ -|- 3a; -f 2 -|- 6a; In a; 
6 (a;- 1)4 


(16) 


With the implementation of the model in the BSM-Toolbox [30] of SARAH [31, 32], we have cross¬ 
checked the one-loop results for both neutrino masses and Br(/r —>■ ey). Moreover, with the SARAH 
FlavorKit [33], we have also checked that the most restrictive lepton flavor violating process in the 
scan to be described below, is just Br(/i —>■ ey). From eq. (9), we obtain 


M 


V 

12 


hla^ah2a constant. 

Q 


(17) 


Comparing this result with the corresponding combination of couplings in the expression for Br(/r —)■ 
ey) in eq. (15), we expect that for a set of hxed input masses Br(/i —)■ ey) turns to be inversely 
proportional to A^. This is illustrated in hgure 2(b) for A = 5 x 10“^, where the scatter plot of 
Br(/i —)■ ey) is shown for the same range of tan/3 values than in hgure 2(a). In such a case, once 
hia are obtained from the Casas-Ibarra parametrization, the specihc hierarchy of Aq £x the several 
contributions to Br(/i —)■ ey). The dispersion of the points is due to the 3 -ct variation of neutrino 
oscillation data [34] used in the numerical implementation of the Casas-Ibarra parametrization, along 
with the random variation of the parameters of 7?.. The minimum value of Br(/i — )■ ey) around 
tan/3 = 1 corresponds to the maximum value of Aq,, while the maximum values happen at the 
cancellation points of each Aq,. In the subsequent analysis, and for a hxed SIM and A, we allow for 
cancellations only by two orders of magnitude from the maximum value of each Aq. 

The full scan of the input masses up to 2TeV, with ms^ > 53 GeV [21] as the dark matter 
candidate, Mr, > 100 GeV to satisfy LEP constraints, and 10“^ < tan/3 < 10^, give to arise the 
dark-gray plus light-gray regions in hgure 3. In particular, the A variation for the SIM with /3 = 7r/6, 
denoted by yellow stars in hgure 2(a), is illustrated with the white dots in hgure 3. The corresponding 
dashed line is obtained for the best-ht values of the neutrino oscillation data and 7?. hxed to the 
identity. The horizontal dotted line in the plot corresponds to the current experimental bound for 
Br(p —)■ ey) < 5.7 x lO”^^ at 90% GL [35]. The upper part of the light-gray region is restricted by 
our imposition to avoid too strong cancellation in Aq . We check that for all the sets of input masses 
in the random scan, this cancellation region always happens when tan /3 < 1. In this way, points 
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Figure 3: Br(/i —)■ 07 ) in terms of the Yukawa coupling strength A for the SIM in eq. (14) with 
(3 = 7 r/ 6 , and the general scan described in the text. 


with tan/3 > 1 are absent from the light-gray region, as labeled in hgure 3. For the same reason, 
in the dark-gray region there are not points with Aq, ^ A,g ~ (a 7 ^ /3 7 ^ 7 ). We can check for 
example that points with Ai ^ A 2 < A 3 are absent inside the dark-gray region of hgure 3. 

The lower part of the dark-gray region is saturated by the values of Mo = 2 TeV, and gives rise 
to the lower bound A > 6 x 10“®. With our restriction in the cancellation of Aq, points in the scan 
with A < 3 X 10“^ can be excluded from the Br(/i — )■ ey) limit. 

5 Collider phenomenology 

The LHC phenomenology in the case of the singlet-doublet fermion dark matter was already analyzed 
in [21]. Their conclusion, is that the recast of the current LHC data is easier to evade, but the 
long-rung prospects are promising, since the region Mp^,mx ^ Mo could be probed up to Mo < 
600 — 700 GeV for the 14-TeV run of the LHC with 3 000 fb“\ 

On the other hand, in the case of the singlet scalar dark matter, the main production processes 
associated with the new fermions remain the same, but there are new signals from the mediation, or 
presence in the hnal decay chains, of the new scalars. The most promising possibility is the dilepton 
plus missing transverse energy signal coming from the production of charged fermions decaying into 
leptons and the lightest scalar. This signal can be important when A is not too large, A < 0.1, and 
3Mn ^ Mo- For a hxed set of input parameters, the random phases in the Casas-Ibarra can be 
chosen to have all the possibilities in the lepton flavor space associated with the coupling hn, with 
i = In view of that, we will focus in the best scenario where Br(x^ —)■ Si) ~ 1. The 

Feynman diagram for the processes is displayed in hgure 4. 

The mass of the charged Dirac fermion can be constrained from dilepton plus missing trans¬ 
verse energy searches at the LHC. In [36], this kind of signals was used by the ATLAS collaboration 
to establish bounds on the slepton masses from the search for pp —)■ —)■ where are 

the sleptons, are the neutralinos and l~ is e~ or p~. Purely left-handed sleptons produced and 
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decaying this way, have been excluded up to masses of about 300 GeV at 95% CL, from the data 
with integrated luminosity of 20.3 fb“^ and the pp collision energy of 8 TeV. This corresponds to an 
excluded cross section of 1.4 fb at NLO calculated with PROSPINO [37]. 

In the present model, the charged fermion held may decay in the mode which are 

proportional to the Yukawa couplings hn. Therefore, a similar hnal state as in the slepton pair 
production is obtained through the process pp —)■ x'^X~ as can be seen in hgure 4. 



Figure 4: Feynman diagram for pp —)■ y+y —)■ . 


In this case, the excluded cross section of this process can be estimated from: 


a{pp /+/ = a{pp x+x ) X Br(x^ ^ (18) 

where a{pp —)■ x^X~) is the pair production cross section of charged Dirac fermion, and Br(x^ —t 
is the branching fraction for x^ —t ^^5'° mode. 

The pair production of charged Dirac fermions can be calculated in the pure-higgsino limit of 
the minimal supersymmetric standard model. The NLO cross section calculated with PROSPINO is 
displayed in hgure 5 as a function of the charged Dirac fermion. 

For points in the parameter space where the Casas-Ibarra solution is chosen such that Br(x^ —t 
Si) ~ 1, and assuming the same efficiency as for the dilepton plus missing transverse energy 
signal coming from left-sleptons in eq. (18), the charged Dirac fermions of the present model can be 
excluded up to 510 GeV, as illustrated in hgure 5. 

Note that many points in the scan of hgure 3 with A < 0.1 and featuring rus^ -C could be 
excluded by this LHG constraint. However, a detailed analysis of the restriction from the Run I of 
the LHG, in the full parameter space of the model, is beyond the scope of this work. 

6 Singlet-doublet fermion dark matter 

In this model, the role of the dark matter particle can be played by either the lightest of the fermions 
Xlop or the lightest of the scalars Si. In the latter case, the present model resembles the singlet scalar 
DM model [1, 2, 3] as long as the other Z 2 -odd particles do not contribute to the total annihilation 
cross section of Si, namely through to the addition of new (co)annihilation channels. Therefore, by 


9 


10^ 



Figure 5: NLO cross section for the charged Dirac fermion pair production at the LHC with pp colli¬ 
sions at a/s = 8 TeV. The horizontal dashed line for the excluded cross section of 1.4 fb, corresponds 
to the mass about 510 GeV illustrated by the vertical dashed line. 


choosing a non degenerate mass spectrum and small Yukawa couplings (which is in agreement with 
neutrino masses) the effects of these particles on dark matter can be neglected. Hence we expect 
that the dark matter phenomenology to be similar to that of the SSDM [38]. 

On the other hand, regarding the case of fermion DM, the present model includes the singlet 
doublet fermion DM model [4, 5, 6, 7, 8, 9]. In such scenario, when the dark matter candidate 
is mainly singlet (doublet) the relic density is in general rather large (small). In particular, a pure 
doublet has the proper relic density for Md ~ 1 TeV [5, 9, 39] with decreasing values as Md decreases. 
Nonetheless, in the present model we have the additional possibility of coannihilations between the 
Z 2 -odd scalars and fermions. In this work, we explore at what extent coannihilation with scalars 
may allow to recover pure-doublet DM regions with M/) < 1 TeV and A < 0.3, while keeping the 
proper relic density. Hereafter, we focus in that specihc region. 

In the simple radiative seesaw model with inert doublet scalar dark matter, the coannihilations 
with singlet fermions can enhance rather than reduce the relic density, as shown in [40]. That work 
also presented a review of the several models [41, 42, 43, 44, 45] where such an enhancement also 
occurs. In particular, supersymmetric models where the neutralino is higgsino-like were considered 
in [45] and it was shown that slepton coannihilations not only lead to an increase in the relic density 
but also to an enhancement in the predicted indirect detection signals. Below, we show that the 
singlet scalars can play the role of the sleptons in our generalization of the higgsino-like dark matter 
with radiative neutrino masses. 

The interactions of the scalars Sa are described by the hia, terms in eq. (1). It turns out 
that Yukawa interactions are suppressed by neutrino masses (hj„ < 10“^) and the same occurs for 
the interaction with the Higgs boson if we impose Af^ < 10“^. In this way the coannihilating scalars 
Sa act as as parasite degrees of freedom at freeze-out leading to an increase of the singlet-doublet 
fermion relic density. 
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Figure 6: Regions consistent with the observed relic density for A = 0.3 and tan/3 = 2. The 
solid cyan line corresponds to the observed relic density without coannihilations which was shown 
to be compatible with the current direct detection bounds from LUX [46] in [9]. The effect of the 
coannihilations with the new scalars is shown for a mass degeneracy of 0.1 to 10% between the scalars 
and the DM candidate. The dark-gray region corresponds to coannihilations with one scalar singlet, 
while the dark plus light-gray regions correspond to coannihilations with two scalar singlets. 


By following the discussion in [40], the maximum enhancement of the relic density is achieved 
when As^ = {ms^ — ''^lop)/’^lop becomes negligible. Accordingly one can write 


^ _ ( 9o + 9s^ y 
go ) 


(19) 


where (D°) denotes the relic density with (without) including Sa coannihilations, gs^ represents 
the total number of internal degrees of freedom related to the scalars participating in the in the 
coannihilation process and go is the total number of internal degrees of freedom when As^ ^ 1. 
When the DM particle is pure doublet (M^j ~ 1 TeV and ^ Tf^)) the fermion masses are 
= Mjv, 3 ~ rn^± = and therefore go = gx 2 + 9x3 + 9x^ ~ Since each real scalar 
have one degree of freedom we have gs^ = 1,2,3 depending on the number of scalars coannihilating 
from which it follows that the maximum enhancement is = 1.27, 1.56, 1.89, respectively. 

This enhancement results in that for the present model with doublet-like DM and A < 0.3 the Md 
required to explain the correct relic density lies in the range [0.9,1.1] TeV instead of taking a single 
value as in the SDFDM model. The values inside this range, arise due to a no mass degeneracy 
between the fermions and scalars. In hgure 6 we show the effect of coannihilations on the relic 
density of m^Qp for a mass degeneracy of 0.1 to 10% between scalar singlets and the DM candidate 
and for A = 0.3 and tan/3 = 2. In particular, in the light-gray region we plot the coannihilations 
with two scalars to facilitate the comparison with the results in [45] for higgsino-like dark matter 
coannihilating with a right-handed stau [g ^ 2 in their plots). As expected, the upper limit in the 
LOP mass is about 20% smaller with respect to the case without coannihilation, and we could expect 

"‘The relic density is calculated with the BSM-Toolbox chain: SPheno 3.3.6 [47]-MicrDMEGAs 4.1.7 [48, 49]. 
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similar enhancements for indirect DM searches as in [45] for g ^ 2. Note that the impact of the 
Sa coannihilations when Mo, Mn < ITeV, is reduced because in such case the dark matter particle 
is a mixture of singlet and doublet (well-tempered DM [50]), and the non-negligible splitting among 
the fermion particles x leads to a non-zero Boltzmann suppression. We have checked that the same 
results are obtained when A < 0.3. 

With regard to DM direct detection in the pure-doublet DM scenario discussed above, it is not 
restricted by the current LUX [46] bounds as long as tan f3 > 0. This is due to the existence of zones, 
known as blind spots, where the spin independent cross section vanishes identically and they occur 
only for positive values of tan (3 [9]^. In consequence, the recovered pure-doublet DM regions are still 
viable in light of the present results of direct searches of dark matter. 

7 Conclusions 

We have combined the singlet-doublet fermion dark matter (SDFDM) and the singlet scalar dark 
matter (SSDM) models into a framework that generates radiative neutrino masses. The required 
lepton number violation only happens if the scalars are real. We have then explored the novel 
features of the final model in flavor physics, collider searches, and dark matter related experiments. 
In the case of SSDM, for example, the singlet-doublet fermion mixing cannot be too small in order 
to be compatible with lepton flavor violating (LFV) observables like Br(/i —)■ ey), while in the case 
of fermion dark matter the LFV constraints are automatically satished. The presence of new decay 
channels for the next to lightest odd particle opens the possibility of new signals at the LHC. In 
particular, when the singlet scalar is the lightest odd-particle and the singlet-like Majorana fermion 
is heavier than the charged Dirac fermion, the production of the later yields dilepton plus missing 
transverse energy signals. For large enough or /i^ branchings, these signals could exclude charged 
Dirac fermion masses of order 500 GeV in the Run I of the LHC. Finally, the effect of coannihilations 
with the scalar singlets was studied in the case of doublet-like fermion dark matter. In that case, it 
is possible to obtain the observed dark matter relic density with lower values of the LOP mass. 
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®Note that tan/3 > 0 corresponds to tan0 < 0 in notation of [9]. 
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A Analytic formulas for masses and mixing matrix 
tral fermions 


The characteristic equation of the mass matrix (2) is [9] 


6 . 


The solutions to the cubic equation in .. are: 


+ 




= 0 


ml =Z2 + 


M, 


N 


=Zi + 


M 


N 


ml =Z3 + 


M. 


N 


where 



q = (MJ + m|) + [MmMI - m|sm(2/3)MB] . 

Notice that g^/4 + p^/27 < 0 and therefore, we have three real masses m^ (i = 1, 2, 
Expanding the eigensystem in eq. (4) by assuming that Nu ^ 0, we have 


where 


= -m?) 

iVii 

+ m, - = -Mh . 

jVli JVij 


Nu = 




- 1/2 


Using the matrix given in the eq. (2), we get the ratios 


N 2 i mx cos 13 Mix — mf) + m\ cos 

Nu mf mf mx{mf sin (3 + Mf) cos(3) 

Nsi [mfiMN — Tnf) + ml cos /3^] 

Nu mx{mf sin/3 + Md cos (3) 

®The analytic formulas for the neutralino masses and the neutralino mixing matrix was analyzed i 


of neu- 


( 20 ) 


( 21 ) 

3). 


( 22 ) 


(23) 

151|. 
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A.l Approximate mixing matrix 

By using the analytical expressions for the mixing ratios of eq. (23) with the approximate eigenval¬ 
ues (5) in eq. (22), we obtain 


Nf, =1 


[Ml + Ml + 2 MdMn sin(2/?)] mj 


(Ml 


Ml? 


+ O (m\) 


^''12 


[sin(2/3) -I- l]m^ 


— — 


2 {M]\j- — Mjy 

[sin(2/?) — l]m| 
2 {Md + Mtv)^ 


+ O (ml) 

+ O (K) 


(24) 


^''21 


m\ (sin f^Mi) + cos 


(Ml 


Ml) 


+ O (ml) 


1 m|(sin f3 + cos (5) [cos /3Mjv — sin {3 {Mn — 2Md)] 

2 AMd {Mn - Mof 

1 ^ m^(cos/9 — sin/3) [sin/? (2 M£)-|-Mat)- f-cos/^M tv] 

2 ^ AMd {Md + Mjv)^ 


+ O (ml) 

+ O (ml) . 


(25) 


^''31 ~ 


Md cos Y -\- Mjv sin (3 

Ml-Ml 


ml + O {ml) 


_1 [MArsin/3 - (Mat - 2M£,) cos/?] (cos/3 + sin/3) _2 , /o 

^ 2 AMd{Mn-Md) ^ 


—- 

^''33 


[Mn sin /3 -|- {Mn + 2Md) cos /3] (cos (3 — sin /3) 
AMd {Mn + Md? 


ml + O {ml) . 


(26) 


In particular, with eq. (5) and the expressions for iV|-, the identity (7) is satished up to terms of 
order O (m^). 
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